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Abstract 

We construct a gravity dual of a boost-invariant flow of an jV" = 4 SU{N) supersym- 
metric Yang-Mills gauge theory plasma with chemical potential. We present both a first- 
order corrected late-time solution in Eddington-Finkelstein coordinates and a zeroth-order 
solution in parametric form in Fefferman- Graham coordinates. The resulting background 
takes the form of a time-dependent AdS Reissner-Nordstrom-type black hole whose hori- 
zons move into the bulk of the AdS space. The solution correctly reproduces the energy 
and charge density as well as the viscosity of the plasma previously computed in the 
literature. 
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1 Introduction 

In the recent years, the application of the AdS / CFT correspondence [1] to the quark-gluon 
plasma (QGP) has become a very active research area. One line of research within such 
holographic studies was initiated by Janik and Peschanski [2] who established a time- 
dependent gravity dual of the boost-invariant flow of an A/" = 4 plasma. This geometry 
has mainly been studied in the regime of large proper time, when the system is near 
equilibrium and approaches the hydrodynamic regime (see however [3l H]). In [5] |12] 
higher-order corrections to this late-time background were constructed and found to be 
equivalent to a gradient expansion of hydrodynamics, see [13] for a review. 

An important aspect of the plasma which has not yet received much attention in 
a time-dependent gravity background is the effects of chemical potentials, even though 
an asymptotic boost-invariant geometry [without corrections) dual to an A/" = 4 plasma 
with f/(l) i?-charge is known for quite some time [H]. Also the transport coefficients 
of plasmas with U{1) currents have already been holographically computed in [T5] |19] 
(up to second order). Such currents are generated, for instance, shortly after the collision 
of two heavy ions, when the two sheets of color glass condensates have passed through 
each other and longitudinal color electric and magnetic flux tubes are produced between 
the sheets [20]. This gives rise to a large topological charge density Fj^'^F^^, which in 
turn leads to an imbalance of the number of quarks with left- and right-handed chirality 
and chemical potentials and /i/,. In addition to the usual baryon chemical potential 
/i = ^^"^^^ , one may therefore also consider a chiral chemical potential /is = which 
mimics the effect of an imbalanced chirality. 
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In this paper we will construct a modification of the Janik-Peschanski background, 
which will additionally include a time-dependent U{1) gauge field. The bulk theory will 
be five-dimensional Einstein-Maxwell gravity with a negative cosmological constant and a 
Chern-Simons term. As in the case without chemical potential, it appears to be difficult 
to find an analytic solution for all times and we will restrict to solving the equations of 
motion at late times. As a further simplification, we seek for a solution in which only the 
time-component of the U{1) gauge field dual to the chemical potential is non-vanishing 
(the spatial components are set to zero). Asymptotically, at large proper time r, we may 
expand the late-time geometry in powers of r"^/^. Employing both Eddington-Finkelstein 
and Fefferman-Graham coordinates we present the late-time solution up to first order (in 
The resulting background will essentially take the form of a time-dependent AdS^ 
Reissner-Nordstrom solution whose inner and outer horizon move into the bulk of the 
AdS space. This background can be extended to a full type IIB supergravity solution 
(by taking the product with an S^) and is dual to a strongly-coupled J\f = 4 SU{N) 
supersymmetric- Yang-Mills plasma with a non-vanishing chemical potential. 



2 Late-time background in Eddington-Finkelstein co- 
ordinates 

In this section we are interested in finding a late-time gravity dual of an expanding A/" = 4 
viscous plasma with non-vanishing chemical potential. 

The relevant five-dimensional Einstein-Maxwell-Chern-Simons action is given by 

S = j d'x (^R+12- F^pF^^ + ^e"^^^'A„F^pF,^ , (2.1) 

where a, /3, ... denote the 5D bulk coordinates. The cosmological constant is A = —6 and 
the Chern-Simons parameter is fixed as k = — l/(2v^). Also, l/(167rG'5) = NI/{%tt'^) for 
an TV = 4 plasma [15]. The corresponding equations of motion are given by the combined 
system of Einstein-Maxwell equations. 

Rap - - Qgap + 2 (^Fa^F^i, - ^Qa^F^^ = , (2.2) 

and covariant Maxwell equations (with Chern-Simons-term), 

V^jF^^" + Ke""^^^" Fp^Fs^ = . (2.3) 
is the field strength of the U{1) gauge field we wish to introduce in the background. 

2.1 Boosted black brane solution 

Our starting point for the construction of a time-dependent solution is the static AdS^ 
Reissner-Nordstrom (RN) black- hole solution [2T]. Using ingoing Eddington-Finkelstein 
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coordinates, we may write the RN metric and gauge field as 
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ds^ = -r2(l - — + ^)dv'^ + 2dvdr + r'^dx^ , (2.4) 

A=-^df, (2.5) 

with mass m and charge q. Here u is a time-hke coordinate (not to be mixed up with the 
scahng variable v introduced below), x are the spatial coordinates on the boundary, and r 
parameterizes the holographic direction. The location of the outer horizon r+ = r+(m, q) 
is given by the largest real positive root of V{rj^) = — mr^ + = 0. 

A charged black hole is dual to a fluid at finite temperature T and chemical potential yU. 
Both the Hawking temperature and the chemical potential are given in terms of r+ by [21] 



These relations can be inverted to give m and q as functions of T and fi |15j . 



m = r^-^, q = ^r^, (2.7) 



with 



Following [HI [To], we now consider the corresponding five- dimensional boosted charged 
black brane solution given by 

1 J + \ ) u^u^dxi'dx" + r'^P^ydx^dx" , 

A = ^^^u^dx^ , Pf,u = + u^Uy , (2.9) 

where is the boost velocity along (/i = 0, 1, 2, 3), and m = m{fi, T) and q = q{fi, T) as 
given by (12.71) . From this solution we may deduce a time-dependent solution by choosing 
the frame = (1,0,0,0) and introducing an Eddington-Finkelstein proper time-like 
coordinate f and rapidity-like coordinate y. We also substitute the asymptotic late-time 
behaviour of T [22] and fi, 

T = Af and /i = /tof"^/^ , A, /io = const. , (2.10) 

into the explicit expressions for m and q. Here we assumed yU oc T, as one would expect 
for a perfect fluid, such that the quotient /i/T = flo/A = const, is independent of time. 
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This leads to the following metricll] 



da^ = -r-Hl - ^ + ^)df2 + 2dfdr + (1 + rffdy' + r^'dxi , 



A=-^df, (2.11) 
with coefficients 

m{f) ^ b(f)-* ^ r^fr?!^ , ,(f) = -^r+(f)^ (2.12) 

For q = (or /io = 0), this metric reduces to the uncharged (zeroth-order) late-time 
solution in Eddington-Finkelstein coordinates found in [ini [HI [12] {'m = = vr^A^f"^/^ 
there). Note that the size of the outer (and inner) horizon r+ (r_) decreases with time. 



2.2 Zeroth-order solution and first-order correction 

The boosted metric (12. lip is not an exact solution of the Einstein-Maxwell equations. It 
is a good approximation of the boost-invariant solution at large f though. At smaller f, it 
receives subleading corrections corresponding to higher-order gradient corrections to the 
energy-momentum tensor and U{1) current, which will be discussed in section [231 These 
corrections to the metric (12. lip can be found by choosing the following metric ansatz for 
the time- dependent solutiono 

ds^ = -r2e"(^'^)c/f2 + 2dfdr + (1 + rffe''^^''^'^dy^ + r^e^^^'^'^dxl , 
A = d{f,r)df. (2.14) 

As in the case without chemical potential, we may introduce the scaling variable v = rf^/'^ 
and expand the metric coefficients in powers of f~^/'^, 

ga{f,r) _ ^^^^ ^ ai(^;)f-2/3 + . . . , (2.15) 
g6(f,r) _ exp(6i(t;)f^2/3 ^ _ _ ^ ^2.16) 
gc(f,r) _ (ji^^^ exp(ci(i;)f-2/3 + (2.17) 

Similarly, for the coefficient of the gauge field we choose 

d(f , r) = D{v)f~^''^ exp(rfi(t;)f + ...). (2.18) 

Note that the gauge field has an overall factor f~^/^. The existence of a late-time scaling 
variable v will be shown in section 13.21 

^ There is an additional 1 in the factor (1 + r-r)^ in front of dy^ which is not expected from the boosted 
sohition (|2.1ip . This is to ensure an asymptotic AdS space in the hmit A — ^ 0, see [TT] for details. 

^For this particular ansatz, the Maxwell equation reduces to -^^^di3{^J~gF^'^) = 0. The Chern-Simons 
term is absent, since only Frf and Ffr are no n- vanishing. 
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The system of Einstein-Maxwell equations f l2.2p and Maxwell equations fl2.3p can then 
be solved order by order in f~^/^. At zeroth-order in f~^/^, we find the coefficients 

= + BW = CM = 1, BW = -^, (2.19) 

where we defined the time-independent variables 

mo = 60 ' = f^^Mr) , go = ) (2.20) 

with m(f) and q{f) as in fl2.12p . In the same way, we also define the variable 

C+ = r+(f)fi/=^ (2.21) 

from the (outer) horizon r_|_ as given by fl2.13p . A{v), ...,D{v) are in agreement with the 
metric (12. lip deduced from the boosted black brane. 

At first order in f~^/^, we find the coefficients 

4g2 2mo C2 

Sv' Sv'^ 

^^r^^^V" 3Cilog[^-0] -41og[t;-0]C' ^9 99^ 

bi{v) = -2ci{v) = -— + C3 + -> —4 , (2.22) 

3v 6 ^ 3Cf - mo 

di{v) = -^ + lv^C^, 
3v 2 

where Q are the solutions of 

Cf-moCf + go' = 0. (2.23) 

The resulting expression for bi{v) is real, even though we need to consider all six roots 
of (I2.23P including the imaginary ones. Explicit expressions for these roots can be found 
in appendix \M Note that one of the six roots of this equation corresponds to the outer 
horizon In Reissner- Nordstrom solutions there is always an upper bound on the 
charge %, at which the discriminant of the equation (I2.23P vanishes, 



% < Qo''- = \l^<- (2-24) 

For larger values of %, there would be a naked singularity at the origin. Remarkably, this 
bound is satisfied for any value of the quotient /io/A and saturated in the limit /io/A — )■ 00, 
as can be seen by substituting (I2.20p with (I2.12p into the bound (I2.24p . In other words, 
there is no bound on the chemical potential. Nevertheless, let us assume that /io <C A in 
order to avoid potential stability problems [23], which arise when the black hole is close 
to extremality. 
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We still need to fix the integration constants Ci^2,3- Ci can be found by requiring 
regularity of the first-order solution fl2.22p at the outer horizon, i.e. Ci should be a function 
of the positive root More precisely, by choosing 

Ci = Id , (2.25) 

we cancel the log[f — C+] terms in bi{v), which are singular aX v = (+. The metric then 
still contains singularities but they are hidden behind the outer horizon. 

The constant C3 is fixed by the requirement that the metric reduces to a pure AdS 
space in the limit A — )■ 0. This simply sets C3 to zero, 

C3 = . (2.26) 

There is one remaining integration constant C2 which can not be fixed at first order. 
Note that, in general, at each order k there is one integration constant which can only 
be fixed by regularity at order k + 1 [10], C2 in our case. Nevertheless, we may guess 
the correct value for C2 by comparing with the uncharged solution flUi [TTl [12] , in which 
C2 = |C|f = I^T'^A^. As for Ci, it seems natural to replace the horizon (h of the uncharged 
solution by the outer horizon of the charged solution such that 

C2 = la . (2.27) 

Later in section [2^3] we will justify this value again. It will turn out to correctly reproduce 
the expected transport coefficients. 

We have checked that for go = (or, equivalently, flo = 0) the metric reduces to 
the first-order corrected uncharged solution found in [101 [HI [12]. Moreover, for the 
Kretschmann scalar we find 

2 4(127gt - 90moqy + ISm^ + 2ggi;6 ^ 10^12) 

^^^upa - :^2 (^-^^^ 



8(254g^ - 150moq^v^ + 24:mlv^ + 2qiv^ - 4:5q^v^C2 + 18moV^C2) _2/3 



^13 



+ ... 



which is only singular at v = 0. In the limit go ~^ 0, we have mo — ?■ vr^A^ and B?^yp„ 
reduces to the corresponding expression in the uncharged case, see |llj . 

We have thus constructed a natural extension of the first-order corrected boost- 
invariant plasma geometry of [lOl [11], [12] to the corresponding one with non-trivial U{1) 
gauge field. 

2.3 Transport coefficients from the background 

In the hydrodynamic approximation, the energy-momentum and U{1) current are given 
by 

{T^^) = ^{Au^u^ + r]^^) + Ii^y, (J^,) = pu^ + T^, (2.29) 
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where the first terms on the right hand side correspond to a perfect fluid with chemical 
potential. Since the velocity field u^, energy density e and charge density p vary slowly 
with the spacetime coordinates, the energy-momentum tensor and current receive higher- 
order gradient corrections given by (up to first order) 

n^. = -r/a^. , = + ^e/'^^u.d^u^ , (2.30) 

where rj, cr and ^ denote the viscosity, conductivity and vorticity coefficient, respectively. 
The corrections satisfy u'^T^, = and u'^H.^v = 0. The transport coefficients of the fluid 
entering these corrections were holographically computed in [151 CBl [I^ (up to second 
order) by slowly varying m'^, q and m in the boosted solution (12.91) with the space-time 
coordinates x'^. In this way the hydrodynamic equations are obtained from AdS/CFT 
without constructing an explicit solution. 

In the following we will compute the first-order corrections directly from our time- 
dependent solution using holographic renormalization techniques [Mj. Recently, a rigorous 
holographic renormalization of the Einstein-Maxwell-Chern-Simons theory, including the 
full back-reaction of the gauge field, has been performed in [25]. The energy- momentum 
tensor can be obtained from 



(T^u) = lim 



f 1 



(2.31) 



where 7^1, is the induce metric on a constant-r hypersurface, which regularizes the bound- 
ary, i^^j, is the extrinsic curvature of on this hypersurface, K the corresponding scalar 
K = K^y'^'^^ and Gf^i, the boundary Einstein tensor with respect to the metric 7^1,. Sub- 
stituting our explicit first-order solution into f l2.3ip . we find the time-dependent energy 
density 



with 



(2^32) 



"° = 8^ ' '^'= 16^^^ - ' ^^-^^^ 



7r2 



and bo = bf~^^^ as in f l2.20p . see appendix [B] for more details on the computation. The 
first term in e{f) is the zeroth-order energy density and is in agreement with that in |15j . 
see Eq. (20a) therein. The second term is the first-order correction and formally agrees 
with that in the uncharged case [lOl [TTl [12] but now with a more general shear viscosity 
^0 = ''7o(a'.0! A)0 This correction is also in exact agreement with the first-order gradient 
correction to the energy-momentum tensor computed in [15]. There [15], the viscosity 
was found to be 



^Asymptotically, f can be identified with the proper time r, f w t, see section |3^ below. 

^In order to check 770 in the limit of vanishing chemical potential, we note that the viscosity is differently 
normalised in [11]. Consider ^ — |C+^o 3^ ^^r jlo = 0. This is identical with rj^'^''^^ = found 
in |llj since ttA — w there. 
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with r+ as in (12 .Sp (The Af = 4 plasma saturates the KSS bound [26]). Here we have 
aheady substituted the asymptotic behaviour T = Af~^/^ and r_|_ = C+f~^/^. Given that 
r/o is defined as ijq = rjf, we get the same rjo as in (I2.33P and thus agreement with [T3] . 
Similarly, the expectation value of the R-charge current can be computed from 

= 4^ (v'-'A'-^^ - ^e^-^-'A^^^Fi^j) , (2.35) 

where Ap"^ is the r"'^ coefficient of the large-r expansion of the gauge field Ap. Since the 
spatial components of the gauge field are zero, the second term proportional to n is absent 
in our case. Substituting the solution for the gauge field into (I2.35p . we read off the U{1) 
charge density 

with % = qf as in dOH]). Recalling ^3^/2 = yur^, we find agreement with the zeroth- 
order charge density in [15], see Eq. (20b) therein. The asymptotic 1/f behaviour of the 
charge density was also found in [T3]. There are no first-order corrections to the charge 
density in our case. 

More generally, for gauge fields with vanishing spatial components, there are no higher- 
order gradient corrections. This follows directly from the relation u'^Ti, = 0. The correc- 
tions T^, are orthogonal to u'^ and cannot come from the near boundary expansion of a 
gauge field proportional to u'^. 



3 Late-time solution in Feffer man- Graham coordi- 
nates 

In this section we seek for a time-dependent solution of the Einstein-Maxwell equations 
(12. 2p and (12. 3p in Fefferman-Graham coordinates. 

3.1 General ansatz and near-boundary behaviour 

In Fefferman-Graham coordinates, we choose the same metric ansatz as in the uncharged 
case [2] given by 

ds^ = 1 (-e"(^'^)rfr2 + e'^^'^^r'dy^ + e'^^'^^dxl + dz') . (3.1) 

Of course, the warp factors a(r, z), b{T, z) and c{t,z) will be modified due to the effects 
from the back-reaction of the gauge field. As before, we set the spatial components of the 
gauge field to zero and assume a non- vanishing time-component, 

Ao = -d{T, z), Ay = A, = A,^ = . (3.2) 
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Let us first study the general behaviour of the solution near the boundary at z = 0. 
Following [3], we choose the small- 2; expansions 

a(r, z) = -£:(r)z^ + ae{T)z^ + as{T)z^ + ... , 

b{T, Z) = h,{T)z* + h,{T)z^ + h^{T)z^ + ... , 

c(r, z) = c^{t)z^ + C6(r)/ + c^{t)z^ + ... (3.3) 

and 

d(r, = p{t)z^ + d4(r)z' + ^(r)^' + ... (3.4) 

Here the lowest coefficients are determined by the energy and charge density, respectively. 
For instance, solving the Einstein-Maxwell equations to lowest order in we obtain 

64(r) = - [sir) + Te\r)) , c^{t) = ^(r) + ^rs'ir) , (3.5) 

as in [5]. There is no back- reaction of the gauge field on the geometry at this order 
(p(r) does neither appear in 64(r) nor C4(r)). Likewise, the metric does not enter the 
Maxwell equations at this order. However, other than the energy density e(r), which can 
be freely chosen (at least at early times), the charge density p(r) is uniquely fixed by the 
2;-component of the Maxwell equations, 

-Ml)_2p'(r)=0, (3.6) 



which is solved by 



p(r) = , (3.7) 
r 



go = const. Any dependence on the warp factors has dropped out in the Maxwell equations 
such that p(r) is independent of ^(t). The result (13.71) for the charge density holds for all 
times r > 0. Remarkably, the charge density diverges at r = 00 

Solving the system of equations (12. 2 p and (12. 3 p order by order, we find the solution 
up to order z^, 



6 ' ' 6 16 ' 128r3 128r2 

' :e^ Hr) + ^^^-^ I z° + 



64r 384 ' ' 36r2 
, e{r), , (2p{rf e'{r) e"{r) 



d{r^ z) = p{r) {z^-^z^+{^-^-^]z^ + ...] {3.1 



^Generic solutions of viscous fluid dynamics are not expected to be regular in the infinite past (see 
footnote 4 in |15j in this context): The volume element on the boundary at constant proper time scales 
linearly with t. Integrating the charge density (oc l/r) over this volume element (cx r) yields a constant 
total charge. Thus, even though the charge density is divergent, the total charge is regular, even at r = 0, 
ensuring the validity of the hydrodynamic approximation. 
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and similar expressions for b{T,z) and c{t,z). These expressions for the warp factors 
generahse the corresponding ones for qo = found in [3] . They describe the all-time near 
boundary behaviour of the background as a function of the energy and charge density. 



3.2 Late-time ansatz for the background 

A full analytical all-time solution is difficult to find, even in the uncharged case (go = 0). 
It is however possible to find a late-time solution. The general late-time behaviour of the 
energy and charge densities can be found as follows (For the energy density the derivation 
is very similar to that in [21 [5]). In the local rest frame the energy-momentum tensor is 
diagonal with elements T^-, Tyy and Txx = ^X2X2 = T^^xs and the current has only a time- 
component Jr while Jy = = Jx3 = 0. Moreover, we assume that these components 
depend only on r. 

Using proper time and rapidity coordinates in fiat Minkowski spacetime, defined by 
Xq = Tcoshy and Xi = rsinhy, 

ds^ = -dr'^ + T^dy'^ + dx\, (3.9) 

the tracelessness condition T^'j, = 0, energy-momentum conservation T^'^ .^y = and charge 
conservation J".,^ = have the form 

-Trr + ^Tyy + Txx = , (3.10) 

rdrTrr + + ^Tyy = , (3.11) 

TdrJr + Jr = 0. (3.12) 

Here we assumed that the anomaly in the f/(l) current is absent, which is true for our 
simple ansatz of the gauge field. 

Comparing with the zeroth-order energy-momentum tensor and current given in f l2.29p . 
in the frame u'^ = (1, 0, 0, 0) we obtain 

<^) = ^^ M^) = 7- (3-13) 

We observe that the asymptotic charge density f l3.13p is in exact agreement with the 
expression (13. 7p for the charge density, which is valid for all times. In other words, the 
late time charge density f l3.13p does not receive any higher-order gradient corrections, in 
agreement with our findings in the previous section. 

Substituting the asymptotic behaviour f l3.13p into the general solution (13. 8p and ex- 
panding the resulting expressions for large r, we get (^o = 1) 

z' 2 + 30goV^/3 6 ^ 10 - 27goV^/3 - 54t'/' 

Q(^^^) = -7i7I+ 27rW3 ' + ' ^ ^ ^ 

z^ 14 + 18r^/3gg -130 + 243r^/3gg-162r«/3 
''l^' - 3^4/3 81rio/3 2916ri6/3 ^+---, (^-15) 

c(t, z) = —j:r + { + ^^TTTT^ Z^ + ■ ■ ■ , 3.16 

^ ' ' 3r4/3 \81rio/3 9^2 / 2916ri6/3 ' ^ ' 
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*..)^^.^-37^.«+(^.|^) + (3.7) 

We find that the dominant terms at large r scale as 

a^{r)z-^-^^^ dn{r)z---^-^, (3.18) 

and similarly bn{T)z"' and c.„(r)2". As in [2], it is therefore useful to introduce the scaling 
variabl^ 

^ = -^. (3.19) 



r 



This suggests the following ansatz at late times, 

a(r, z) = ao{v) + ai{v)^^ + ... 

rf(r, z) = (^d,iv) + rfi(i;);^ + (3.20) 

and similarly for b{T,z) and c(r, 2). Inserting the ansatz (13.11) and (13. 2p with (I3.20p into 
the combined system of Einstein-Maxwell and covariant Maxwell equations (12 .2^ and (12. 3 p 
will turn the equation of motions into a system of nonlinear ordinary differential equations 
for the coefficients ai, ...,di {i > 0). In principle, this system can then be solved order by 
order in r~^/^. 



3.3 Zeroth-order solution 

In the following we restrict to give an exact solution for the zeroth-order coefficients 
ao{v), ...,do{v). The non- vanishing components of the Einstein-Maxwell equations are 

(rr) : Ae-^^^'^h^^do'iv)^ = Qbo'{v) - vh^{vf + 12co'(tO 

-2vho{v)co'{v) - 2vco'{vf - 2vbo"{v) - 4:VCo"{v) , 

(yy) : 4:e-'''>'^''h^do'{vY = -Gao'iv) + vao'{vy - 12co'{v) 

+2vao'{v)co'{v) + 3vco'{v)^ + 2vao"{v) + 4:VCo"{v) , 

(±±) : 4e-'^«(^)t;3do'(t;)2 = -6ao'{v) + vao'ivf - 6bo{v) + vao'{v)bo'iv) + vboivf 

-6co'{v) + vao'{v)co{v) + vbo'{v)co{v) + vco'iv)'^ 
+2vao"{v) + 2vbo"iv) + 2vco"{v) , 

{zz) : 4e-'^«(^)t;3rfo'(f )2 = 6a'o{v) + Gb'^iv) - va'f^{v)b'Q{v) 

+12c'oiv) - 2va',{v)c',{v) - 2vb',iv)c'^iv) - vc'^ivf , 

(zt) : Qa',iv) - Ab'.iv) - va',{v)b',{v) + vb',{vf + Ad,{v) 

-2va'Q{v)c'Q{v) + 2vc'Q{vf + 2vb'l^{v) + Avcl{v) = . 

(3.21) 

^With hindsight, this justifies the introduction of the scaling variable v = rf^^^ in the previous section 
for the late-time solution in Eddington-Finkelstein coordinates. 
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At zeroth order, the z- and r-components of the Maxwell equation both lead to the same 
equation, 

(2 + mo'(f ) - vh^{v) - 2vc^{v)) do\v) = 2vdo"{v) . (3.22) 

The other components are zero. 

These equations can be simplified a lot. Note that only four out of the five plus one 
equations are independent. We also find from a linear combination of the rr- , zz- and 
2;r-components of the Einstein-Maxwell equations that bo{v) = Cq{v). Next, the Maxwell 
equation fl3.22p can be solved for do{v), 

V 

do{v) = Si f ve^Mv)-bo{v)-2co{v)) ^-^ ^3^23) 



where 5*4 is some integration constant which will be fixed below. 

Substituting this back into the Einstein equations, the two remaining independent 
equations are given by the rr- and zz-components. The first one (rr) is an equation for 

SvKiv) + 3v{b'Q{v) f - %{v) + 8g2|;5e-3boM = g . (3.24) 
while the second one {zz), 

can be used to find ao{v) as soon as a solution for bo{v) is known. Our primary goal will 
be to solve 03.241) for bQ^v). ao(i') and do{v) can then easily be obtained from fl3.25p and 

Later, in order to fix some integration constants, we will need the asymptotic solution 
close to the boundary which can be expanded in powers of v as (here we present it up to 
O(i;i0)) 

/ N 4 , IOQ'o^ 6 ^0 8 2go^£o 10 , 
ao{v) = -eov + — - -^v +■■■ , 

. / \ / \ ^04 2go^ 6 ^0 8 I 14go^£^o 10 , /o 

br)(v) = Cn f = — V V -V H V + ■ ■ ■ , [6.Zbj 

uv ; ov ; 3 g ^3^ 

do{v) = qov —V + —V + —^v H . 

It shows us that the solution exists and is uniquely fixed by parameters Eq and go- Com- 
paring the expression (13.231) with the boundary behaviour fl3.26p . we may immediately fix 
the integration constant S4 as 

^4 = 2qo . (3.27) 
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We now solve fl3.24p for &o(^)- By setting 

bo{v) = log{(3{v)) , (3.28) 
we simplify this equation to the form 

v^" - 3/3' + ^qyr^ = , (3.29) 



which turns out to be the modified Emden- Fowler equation [27|. Its solution can be 
written in the parametric form 

(3{v) = pSi exp I 52 I (^^ + ^1 + ^ dp \ (3.30) 



p+ 



and 

p 



- = 53exp<;||(^ + 5, + ||^) dp}. (3.31) 

p+ 

Here 5*1 , S3 and p+ are some integration constants. One can in principle absorb p+ in 
5*3 but we separate them for the moment. There are two useful expressions for f3{v) and 
/3'(^), 

(3{v)=pv^ (3.32) 

and 



From fl3.26p . we get the near-boundary conditions 

^{v) = 1 + 0{v') , (3'{v) = + 0{v') , (3.34) 

which will be used to fix the integration constants 5*1 and 5*2. 

Comparing fl3.32p with fl3.34p . we find that near the boundary v should behave as 

v^^, (3.35) 
which is small if p is large. This should be compared with the general large-p behaviour 

..S3exp^^/(^ + ...V'#H»-*./'. (3.36) 



2 „ 

p+ 
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This fixes So as 



^2 



1 

2 



(3.37) 



Substituting f l3.35p into f l3.34p . we extract the expected asymptotics for as a 

function of p, 

4^0 1 



Generically, at large p, fl3.33p is approximated by 



-4v\ Si 



T 



2pv 



4^1 



p. 



,3/2 



which fixes 5*1 as 



Si 



£o 
' 3 ■ 



(3.38) 



(3.39) 



(3.40) 



Remarkably both constants 5*1 and ^2 do not depend on qq. 

Let us now relate 5*3 and p+ by setting 5*3 = v+ with v+ = v{p+). v+ will be fixed by 
the requirement that the outer horizon of the geometry is located a.t v = Formally, 
the horizon v+ is defined as the largest zero of the denominator on the right hand side of 



(3.41) 



This can be rewritten in terms of p+. Using fl3.33p . fl3.32p and the expressions for 5*1 and 
5*2, we get the condition 

4^0 , 4g2 



p;-— P+ + ^ = U, (3.42) 
which can be solved by Cardano's formula. The largest solution of this equation i^ 

1/3 / / I 7T— \ 1/3 / : l/3^ 



P+ 




(3.43) 



The last step is to fix t>+ in (13.311) . This can be done by substituting the p+ solution 
(13.431) (and all the constants 5*1^2,3) back into fl3.3ip and expand v{p) for large p. In this 
way we determine the constant on the right hand side of fl3.36p as a function of v^. Since 
this constant must be one, we get 



v+{eo,qo) = exp <J - lim 



dp 



logp 



(3.44) 



^ In order to extract the roots correctly we use the fohowing standard convention: 

1/3 



'2^4/3 



£0, where a± = [-q^± J - 
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where U{p) is defined as 



U{p) 



p2 



£o 
3 



3p 



(3.45) 



For go = 0, the integral can be performed analytically and w+ reduces to the well-known 
result for the horizon [2], 



v+{eQ,qo)\ 



(3.46) 




0.2 0.3 0.4 0.5 0.6 0.7 



Figure 1: The outer horizon t>+ = v^^Eq, go)- 

For general go, this integral can in principle be written as a lengthy expression of 
elliptic integrals of the first and third kind, F{(j), k) and n(n; 4>\m), respectively, which we 
will not do here. Instead, in Fig. [1] we show the dependence of w+ on the charge go for 
some particular choices of ^o- We note that for each Eq there is some maximal allowed 
value of the charge at which the black hole becomes extremal. This value can be found 
from the condition that the discriminant A of fl3.42p vanishes. 



A 



which leads to the bound 



qo < Qo 



extr. 



-3/4 



0, 



(3.47) 



(3.48) 



In Fig.[2]we present some plots of the exact solution and compare them with the power 
expansions fl3.26p . For the particular choice Eq = 1 and go = 0.6 < go^**"' the difference 
between both curves is clearly visible. The function ao(w) by definition has a singularity 
on the horizon, as can be seen in Fig 2(a), The other functions 6o(f), c^iy) and d^iy) are 
regular on the horizon and their power expansions are valid up to f < 1. Note also that 
dQ{v) grows quadratically near the boundary, which reflects the Coulomb law in D = 5 
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(a) (b) (c) 

Figure 2: Exact solutions (red curves) and their near-boundary power expansions (blue 
curves) for Eq = 1 and go = 0.6 < qQ^^"^'. The black dashed lines correspond to the horizon 
Vj,{eo,qo) ~ 1.685. 



dimensions. Near the horizon it approaches some finite constant value /xo related to the 
chemical potential as 



- A\ - 4 I - ^0(^+) - ^O(O) _ /^O (O AG] 

H' -^Olboundary ^O\horizon -j-l/S 7-I/3 ' yo.HiVj 



which confirms the scaling behaviour fl2.10p . 

In summary, the zeroth-order solution bo{v) is given by 

CP 



g6o(«) 



V = V. 



1 

pvl exp ( - - / U ipy^^^dp) , 

exp(-i r Uipy'/'dp) , (3.50) 



with f+ and U{p) given by fl3.44p and f l3.45p . respectively. ao(f) and do{v) are obtained 
by substituting bo{v) in fl3.25p and f l3.23p . 



3.4 Feffer man- Graham vs. Eddington-Finkelstein coordinates 

The zeroth-order solution in Fefferman-Graham (FG) coordinates can be related to that 
in Eddington-Finkelstein (EF) coordinates by the coordinate transformation 

f = T, r = -e''("'")/^ (3.51) 
z 

Transforming the Eddington-Finkelstein metric f l2.14p - fl2.18p and comparing the result 
with the power expansion f l3.26p . we find 

go = ^qo- (3.52) 
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Comparing also the bound (13.481) with that in Eddington-Finkelstein coordinates given 
by fl2.24p . we find some relation between Eq and mo using f l3.52p : 



(3.53) 



For qo = this relation can be easily checked by our solution with the general form of the 
metric in [21 [5l [12]. In this case, we find that our solution (I3.50p reduces to 



1 , ^0 4 
1 + -. 



and similarly ao(f), such that 



1 90=0 



mo„,4\2 



(1 - 

1 + ^v'^ 



-dr' 



r'il + '^v') dy'+(l + '^v^]dx', 



(3.54) 



(3.55) 



We also note that the transformation (I3.5ip relates the outer horizons, C+ in EF 
coordinates and p+ in FG coordinates, as C+ = ^/p+. The chemical potential fio can 
therefore be written as a function of Eq and go- Using (I2.12p and (I3.52p . we find 

This dependence is shown in Fig. [3] for some particular values of Eq. 



(3.56) 




0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 



Figure 3: The chemical potential /xo = /^o(^o,?o)- The dotted line corresponds to the 
upper bound for fiQ. 



Using this expression with the definition (I3.42p for and substituting there the 
maximal value for go (13.480 we find the following bound for the chemical potential: 

1/3 



/io(^o,go) < Ato 



extr. 



extr. 



-1/4 



(3.57) 



This is not in contradiction with our earlier statement that the disappearance of the 
horzion does not impose a bound on /io/A. Note that if we identify Eq with tuq as in 
(I3.53p . then Eq explicitly depends on /io = fio and therefore (13.570 is not a bound on /io/A. 
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4 Conclusions 



We constructed a natural extension of the late-time boost-invariant background found 
in [2] (and [5]-[T2]) to a background dual to an expanding A/" = 4 plasma with chemical 
potential. The solution we found depends on two parameters, the chemical potential fto 
and temperature scale A, which are encoded in the mass parameter rriQ and charge go oi a 
time-dependent AdS Reissner-Nordstrom-like solution. In Eddington-Finkelstein coordi- 
nates the first-order solution is given by the expansion fl2.14p - fl2.18p . with the zeroth-order 
and first-order coefficients given by fl2.19p and fl2.22p . respectively. We showed that the 
viscosity of the boundary theory computed from the time-dependent solution is in agree- 
ment with that in [13]. We also constructed a zeroth-order solution in Feffer man- Graham 
coordinates, which we presented in parametric form, see the general ansatz fl3.1l) and fl3.2l) 
with fl3.50p . FG coordinates may be the preferred choice, when strings [2S] or branes 
j29| [30] are embedded into the geometry. Finally, we found the coordinate transformation 
which maps the zeroth-order solution in FG coordinates to that in EF coordinates. 

We argued in several ways that the charge density behaves like at all times. Unlike 
the energy density, it can not be chosen freely at early times. This is basically because 
the charge density behaves like at large r, see e.g. fl2.36p or fl3.13l) . and higher-order 
corrections are absent. It also follows directly from the equations of motion, see (13.70 
which holds for all times. Naive extrapolation to early times shows a singularity in the 
gauge field at r = 0. However, this does not signal a breakdown of the hydrodynamic 
approximation since the total charge is constant at all times and therefore regular even 
at r = (see footnote [S] on p. [TUl) . 

A possible application of the background, when appropriately modified and extended, 
could be the chiral magnetic effect (CME) [31]. The CME states that, in the presence 
of a magnetic field and non-zero chiral chemical potential, an electromagnetic current of 
the type J oc /xsB is generated in the plasma. The CME is a non-equilibrium process 
and requires the introduction of gauge fields with (time-dependent) spatial components. 
For instance, for an electric field E, one needs to introduce the spatial component ^3 = 
As{t, r) with boundary condition — tE at r — > 00. The time-dependence of the gauge 
field reflects non-equilibrium physics and requires the back-reaction on the geometry, 
unless one keeps it infinitesimal [52] (see also [331 [M] for an AdS/CFT approach to the 
CME). In this case one would also obtain higher-order corrections to the charge density 
due the effects of the Chern-Simons term, which are absent in our solution. An attempt 
to include an E field in the dual of an electrified plasma (without chemical potential) has 
been made in [T7] . 

Finally, it would be interesting to find a numerical solution of our background d la 
Chesler and Yaffe [1] which would hold beyond the hydrodynamic regime. 
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Appendix 



A Roots of ( [2:231 ) 



For completeness, we present the six roots of (12.231) in this appendix. The equation f l2.23p 
is depressed bicubic in d and, therefore, can be solved by Cardano's formula. It has six 
solutions, which can be expressed as 



del ±^/a^^+7^, ± 



where 



a 



± 



2 V 4 27 



(A.l) 



(A.2) 



Here we use the standard convention = mo/3. One can recognize the outer horizon 

= r_|_r^/^ in the first pair of solutions and the inner horizon C_ = r_T^^^ in the second 
one. 



B The energy-momentum tensor (12.311) 



In this appendix we introduce the geometric quantities used for the computation of the 
energy-momentum tensor (12.31 



with induced metric 7^1, on itl 



Here we consider an r = const, four-dimensional surface 



(B.l) 



where (7^1, is the 5-metric and is the outward-pointing unit normal vector to the surface. 
For our ansatz (I2.14p it is given by 



n^= 0,0,0,0, 



y/-9ff 



(B.2) 



where gff = —r'^e'^^'^'^\ The indices of the induced metric can be raised and lowered by 
means of the 5-metric gfj_^, 



The surface extrinsic curvature is given by 



(B.3) 



(B.4) 



'Here and after all Greek letters denote a 5-index. 



20 



where we put to covariant derivatives associated with the induced metric, while the 
derivatives on the right-hand side are defined with respect to the 5-metric. We also define 
a scalar K = K^j^^g^'^ = Kfj^^,'-/^'^ , which is used in the Gibbons-Hawking- York part of 

The Einstein tensor on the surface is defined as 

G^. = ^'^R^.-^j^.^^^R, (B.5) 

where the 4-tensors can be expressed through the 5-tensors (defined with respect to (7^1,) 
by the Gauss equations: 

(^)i?^, = (^)i?%;3.7a^= 7/7/7/ i?%A<x + /^i^M.-/^Mai^^, (B.7) 

(^)i? = ^'^^ R^^Y'' = R - n"n^ Raf3 + - Ka/^K'''^ , (B.8) 
where the raising/lowering rule is given by K^^ = K^^l"''^ = luaK°'^. 
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